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We calculate the phase diagram of the SU(iV) Hubbard model describing fermionic alkaline earth 
atoms in a square optical lattice with on-average one atom per site, using a slave-rotor mean-field 
approximation. We find that the chiral spin liquid predicted for N > 5 and large interactions passes 
through a fractionalized state with a spinon Fermi surface as interactions are decreased before 
transitioning to a weakly interacting metal. We also show that by adding an artificial uniform 
magnetic field with flux per plaquette 2i t/N, the chiral spin liquid becomes the ground state for all 
N > 3 at large interactions, persists to weaker interactions, and its spin gap increases, suggesting 
that the spin liquid physics will persist to higher temperatures. We discuss potential methods to 
realize the artificial gauge fields and detect the predicted phases. 


Introduction .—The experimental realization of a topo¬ 
logically ordered phase of matter other than the frac¬ 
tional quantum Hall effect that occurs in two-dimensional 
electron gases is a major goal in both condensed matter 
and atomic physics. Phases with intrinsic topological or¬ 
der pQ are of fundamental interest, as they exist outside 
of the standard symmetry-breaking framework for clas¬ 
sifying phases of matter and display exotic phenomena 
such as fractionalized excitations and edge states that 
are robust to local perturbations [2]; in some cases these 
phases have been predicted to be useful for topological 
quantum computation Ultracold atomic systems 

are uniquely tunable and clean systems that offer a plat¬ 
form to realize exotic phases. However, so far, reach¬ 
ing the required low temperatures remains a challenge. 
Previous work predicted a topologically ordered chiral 
spin liquid (CSL) ground state in fermionic alkaline earth 
atoms (AEA) in a deep square optical lattice [5], 6 . In 
this Letter we show, within a slave-rotor approximation, 
that by applying a synthetic gauge field to this system 
it is possible to enhance the parameter space where the 
CSL exists, to increase the corresponding spin gap, and in 
turn to increase the temperatures at which CSL physics 
manifests. In addition, without a synthetic gauge field, 
away from the strongly insulating limit we find a gapless 
quantum spin liquid with a spinon Fermi surface. 

Recently, experiments have trapped and cooled AEA 
to quantum degeneracy and loaded them in an optical 
lattice tans. Moreover, experiments [T8H20] have con¬ 
firmed the predicted SU(7V) spin symmetry in the colli- 
sional properties of fermionic AEA I2T1I231 . This SU (N) 
symmetry generalizes the usual SU(2) symmetry, and N 
can be controllably varied by initial state preparation up 
to 21 + 1, with I the nuclear spin (as large as N = 10 
for 87 Sr with I = 9/2). The low temperatures reached in 
recent experiments [24H26] , at which short range spin cor¬ 
relations should begin to develop, makes it particularly 
timely to study quantum magnetism in these systems. 


Several theory works have addressed questions related to 
the expected SXJ(N) magnetic phases in the strongly in¬ 
teracting limit mmmm- 

In parallel, other ultracold atom experiments have real¬ 
ized synthetic gauge fields [44H5T] . In these experiments, 
the atoms behave as if they were charged particles in ex¬ 
ternal electromagnetic fields despite their neutrality. Al¬ 
though many schemes in principle can create the gauge 
field that we study in this paper, we focus on methods 
utilizing laser-induced tunneling [52U54] . 

AEA in optical lattices with synthetic gauge fields .— 
AEA in a sufficiently deep optical lattice are described by 
an SXJ(N) generalization of the usual (.N = 2) Hubbard 
model, 

H = ~ t J2 e ^ u 4,i c a,j + \ “ C (!) 

(i,j),OL i 

where c a ^ is the fermionic annihilation operator for nu¬ 
clear spin state a at lattice site i, 12 (ij) indicates a 
sum over nearest neighbors i and j; <f>ij = is 

the (externally imposed) lattice gauge field. We define 
n i — 12 a i c a i> an d ^ and ^ are th e hopping energy 
and on-site interaction energy, whose ratio can be tuned 
by modifying the optical lattice depth. In this Letter, we 
take the average fermion number per site to be one. 

The gauge field <f>ij depends both on the artificial elec¬ 
tromagnetic field as well as the gauge choice. We are 
interested in the physics of a two-dimensional square lat¬ 
tice with a spatially uniform, time-independent artificial 
magnetic field, and use the Landau gauge where 

^Xjd yj -i j y i if {z, j} bond is vertical ^ 
0 otherwise, 

Xj is the x coordinate of site j measured in lattice units, 
and is the flux penetrating a single square plaquette of 
the lattice [55]. We focus on the case <J> = 2i t/N, because 
this choice of $ is favorable for the existence of the chiral 
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FIG. 1. Phase diagram, calculated with a slave-rotor mean-field approximation, as a function of spin degrees of freedom N 
and interaction strength U/t in the (a) absence and (b) presence of an artificial uniform magnetic field with flux per plaquette 
<f> = 2tt/N, illustrated in panel (d) for N = 3. Thin black lines are second order phase transitions, while thick black lines are 
first order phase transitions. The states found are the valence bond solids (VBS), chiral spin liquid (CSL), spinon Fermi surface 
(SFS), Fermi liquid (FL), and integer quantum Hall (IQH) states. These are described in the text and illustrated in panel (c). 


spin liquid. We note that the magnetic unit cell associ¬ 
ated with the translational invariance of the Hamiltonian 
is enlarged from the one imposed by the optical lattice 
potential. Figure [ljd) shows the system with this flux 
and gauge choice, and the enlarged magnetic unit cell, 
for N = 3. 

We calculate the phase diagram and properties of this 
system within a slave rotor mean-field approximation 
[56] [57], which we describe briefly. This technique is de¬ 
signed to match on to the previous large-TV solution in 
the large U/t limit, and is well-suited for describing non¬ 
magnetic ground states in proximity to the Mott tran¬ 
sition. First we expand the Hilbert space to include a 
U(l) bosonic rotor degree of freedom on each site, Oj, 
and new fermionic spinon degrees of freedom associated 
with operators which are defined by 

c a j = e-** f a , j. (3) 

In order to reproduce the original Hilbert space, we must 
impose the constraint 

ij-E < 4 > 

a 

that the rotor angular momentum Lj is uniquely de¬ 
termined by the particle number. Here, Lj satisfies 
[Oj,Lj] = i. We rewrite the Hamiltonian in terms of 


these new degrees of freedom, giving 

H = -tJ2 + I ( 5 ) 

(hj),OL i 

Although the rewritten Hamiltonian Eq. © together 
with the constraint Eq. © is exactly equivalent to 
Eq. ([!]), to make further progress we make a mean- 
field approximation to decouple the rotor and spinon de¬ 
grees of freedom. We then obtain the coupled mean-field 
Hamiltonians for the rotors and the spinons, 

H r = ~J 2 -W' + E + hi[Li + 1} ’ (6) 

(ij) 1 

H f -~ E - E ( 7 ) 

(i,j),a i,a 

where hi is a Lagrange multiplier that enforces on aver¬ 
age the constraint Eq. Q, Lj = t(e l6i ~ l6j ) r , and = 

te 1 ^ J2a(fi,ifa,j)f- Here the sub-index r (/) refers 
to taking the expectation value in the rotor (spinon) 
mean-field ground state \L) r (|V’)/)- The Hamiltoni¬ 
ans H r and Hf are invariant under a U(l) gauge trans- 
formation, /V -> Oi + x». and Uj ->• 

tije lXi ~ lx L Jij Jije~ lXi+lXj . We solve H r and Hf self- 
consistently for several variational ansatz [5^ and find 
the ground state by optimizing the total energy ( , ijj\H\'ijj) 
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where H is given by Eq. © and I'lp) = |'0) r |V ; )/ is the 
mean-field state. 

Results .— Figure [lja, b) shows the slave-rotor mean- 
field phase diagram as a function of U/t and TV; the top 
panel shows the phase diagram in the absence of a gauge 
field and the bottom shows the phase diagram for a gauge 
field with flux = 2tt/N. We find five phases: Fermi 
liquid (FL), integer quantum Hall (IQH), valence bond 
solids (VBS), a gapless spin liquid with a spinon Fermi 
surface (SFS) [59]. and a chiral spin liquid (CSL) 6(11 ITJTj . 
Thin black lines indicate second order transitions and 
thick black lines indicate first order phase transitions. 
Generically, the role of the Hubbard U interaction is to 
localize the atom on lattice sites. Such Mott localization 
is signalled in the rotor sector; when the bosonic rotor is 
gapped and uncondensed with (e l9 ) = 0, the system is in 
a Mott insulating state. The mean-field parameters and 
some key properties of the different phases are listed in 
Table [TJ As we show in the table, the rotor and the spinon 
may experience different, even opposite, gauge fluxes in 
their mean-field Hamiltonians for different phases. Since 
the rotor and the spinon must form a whole atom, the 
total gauge flux experienced by the rotor and the spinon 
should be equal to the synthetic gauge flux that is exter¬ 
nally imposed on the atom. 

The FL phase is very similar to the usual SU(2) Fermi 
liquid, and its structure and instabilities are essentially 
those described in the absence of a lattice |22]. The VBS 
are translation-symmetry breaking phases with repeating 
units of SU(TV) singlets spread across multiple sites. In 
particular, as we plot in Figure [ljc), the system is decou¬ 
pled into 6-site rectangular (4-site square) clusters in the 
SU(3)-VBS [SU(4)-VBS] state. The SFS spin liquid state 
is characterized by a gapless spinon Fermi surface with a 
gapped bosonic rotor in the mean-field theory. Going be¬ 
yond the mean-field description, we need to include the 
U(l) phase fluctuation of the spinon hopping tij. This is 
the internal gauge fluctuation E3; it is dynamically gen¬ 
erated and is unrelated to the synthetic gauge field that 
is imposed externally. At low energies, the SFS spin liq¬ 
uid is described by the spinon Fermi surface coupled by a 
fluctuating internal U(l) gauge field [57l I62H661 . Due to 
the spinon-gauge coupling, the overdamped U(l) gauge 
fluctuation scatters the spinons on the Fermi surface and 
destroys the coherence of the spinon quasi-particles. The 
resulting state is a non-Fermi liquid of fermionic spinons. 
The CSL is distinct from the SFS in that the spinons 
form an integer quantum Hall state in the CSL. Upon 
coupling to U(l) gauge fluctuations, this leads to a chiral 
topologically ordered phase with anyon excitations, and 
gapless chiral edge states that carry spin but no charge 

m • 

To understand the global structure of the phase dia¬ 
gram, it is useful to consider the two limits U/t = 0 and 
U/t oo. The FL and IQH states are simply the non¬ 
interacting ground states occurring at U/t = 0. In the 


Phases 

(e i$ ) 

rotor flux 

spinon gap 

spinon flux 

FL 

/0 

0 

0 

0 

SFS 

0 

0 

0 

0 

CSL 

0 

-2tt/TV 

^0 

2tt/TV 

SU(3)-VBS 

0 

—7r 

^0 

7 T 

SU(4)-VBS 

0 

0 

^0 

0 

IQH 

/0 

0 

^0 

27 t/N 

CSL 

0 

0 

^0 

2t t/N 

SU(3)-VBS 

0 

7 t /3 

^0 

7r 

SU(4)-VBS 

0 

7 t /2 

^0 

0 


TABLE I. Parameters that characterize the obtained phases. 
The upper five (lower four) rows describe phases in the ab¬ 
sence (presence) of the synthetic gauge field. The rotor 
(spinon) flux refers to the flux that is experienced by the rotor 
(spinon) in the mean-field Hamiltonian H r ( Hf ). For the FL, 
SFS, IQH, and CSL states, the flux is defined for the elemen¬ 
tary square plaquette. For SU(3)-VBS [SU(4)-VBS] state, the 
flux is defined through the 6-site [4-site] cluster [58] . 

strongly interacting limit, the Hubbard model reduces 
to an SU(TV) Heisenberg model, and the phase diagram 
coincides with previous slave-fermion mean-field calcu¬ 
lations of the Heisenberg model [5]: for TV = 3,4 the 
ground state is a VBS, while for TV > 5 the ground state 
is a CSL. This is true both with and without a synthetic 
gauge field, as in the the U/t —>> oo limit the physics 
is governed by two-site nearest neighbor superexchange, 
which is insensitive to the gauge flux. 

In the intermediate U/t regime, the gauge field causes 
more significant differences. Without a gauge field, we 
find that an SFS phase intervenes between the non¬ 
interacting FL and Heisenberg-limit CSL or VBS for all 
TV except TV = 4, in which case there is a direct transi¬ 
tion between the FL and VBS ground states. The FL- 
SFS transition is second order and is expected to remain 
continuous beyond mean-field theory [67], while the SFS- 
CSL and FL-VBS are first order phase transitions. In 
contrast, in the presence of the 4> = 2tt/N gauge flux, 
a direct second order transition occurs between the non¬ 
interacting IQH phase and the CSL phase within our 
mean-field theory, and the CSL exists at intermediate 
U/t even for TV = 3 and 4. 

The gauge field increases the parameter space for which 
the CSL occurs: in addition to persisting down to TV = 
3,4, the CSL occurs for a broader range of U/t values. 
In particular, the minimum U/t for which the CSL exists 
decreases from about U/t « 5.5 to U/t « 3.5 (the exact 
values depend on TV). 

In the CSL, both the spinon sector and the rotor sector 
are gapped. Figure [2] illustrates the excitation gap A’s 
dependence on U/t, TV, and the gauge flux in the CSL 
where A is the smaller of the spin gap and the rotor 
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FIG. 2. The excitation gap of the CSL phase, A, as a function 
of interaction strength, U, both in units of the tunnelling t . 
The curves illustrate the N- and magnetic flux ^-dependence. 
From bottom to top, we show (N — 10, <F = 0); (N = 10, <f> — 
2 tt/10); (N = 5,<f> = 2 tt/ 5); and (N = 5, <f> = 2 tt/ 5). The 
turning points at U/t w 4 are the locations below which the 
rotor gap becomes smaller than the spin gap. 

gap. In the slave-rotor mean-field approximation, the 
spin gap is simply the band gap of the spinon spectrum, 
and the rotor gap is set by the Hubbard U interaction 
and thus stays much larger than the spin gap in the Mott 
insulating regime except near the Mott transition. For 
a given U/t, the spin gap slightly increases when the 
gauge field is turned on. An even more favorable effect 
of the gauge field for the spin gap occurs because the 
CSL persists to lower U/t. Since A increases as U/t 
decreases, the gauge field increases the maximum A by 
about a factor of 1.5. Because A sets the temperature 
to which the CSL’s characteristics remain, we therefore 
expect the gauge field to increase the temperature range 
over which the CSL behavior is accessible. 

Gauge field implementation. —Many proposals to im¬ 
plement artificial gauge fields exist. Here we suggest one 
scheme, which uses Raman-induced tunneling in deep lat¬ 
tices subject to a uniform potential gradient HU[50]. A 
Raman process is on resonant with the energy splitting 
between adjacent lattice sites, and the atoms acquire a 
phase kick each time they hop, imprinting the phase faj 
in Eq. 0 - This scheme is natural for our current consid¬ 
erations, since it utilizes the optical lattice and generates 
the Hamiltonian Eq. 0 with strong gauge fluxes. Gauge 
fields have been recently demonstrated in bosonic alkali 
atoms using this technique 49. [50], although we note 
that these experiments have observed unexplained heat¬ 
ing, which could be problematic for realizing low temper¬ 
ature phases. 

We also mention the alternative scheme proposed in 
Ref. [54] that seems natural for the present work with 
AEA: rather than using Raman lasers, one traps the 1 So 
ground (g) and 3 Po excited (e) states in, for example, 
a checkboard pattern in an optical lattice by using an 
appropriate, “anti-magic,” wavelength [68]. The e state 


has a ^ 100s natural lifetime, and is therefore stable on 
the timescale of the system. Because a single laser can 
directly drive tunneling of a g atom to an e atom at an 
adjacent lattice site while imprinting a phase faj , one 
avoids the complexity of driving Raman processes. How¬ 
ever, when this proposal is implemented in the context 
of interacting quantum phases additional considerations 
arise that were not accounted for in the prior analysis. 
First, two e-state atoms on the same site can inelasti- 
cally collide and be lost from the trap. We have found 
that this problem can be largely mitigated when using a 
checkerboard g-e pattern [69] . Second, the interactions 
are inhomogeneous, being different for the sites occupied 
by g atoms and e atoms. This issue can modify the dis¬ 
cussed phase diagram. Third, the flux generated in the 
simplest implementation of this proposal is staggered and 
thus requires rectification techniques to make it homoge¬ 
neous. 

Preparation and detection. —Reaching the temperature 
regimes to observe the phase diagram Figure [l] is chal¬ 
lenging. However, the expected advantage of the SU(7V) 
symmetry for cooling PI H3 ES CD] together with the 
less stringent temperature requirements to observe CSL 
phases in the presence of the synthetic gauge field might 
help achieve the required conditions. Other potentially 
favorable aspects of the gauge field are the absence of 
an intermediate SFS phase and that all transitions are 
second order in the mean-field analysis. Consequently, 
adiabatically going from weak to strong interactions may 
be easier than in the absence of the gauge field. On the 
other hand, the gauge field itself introduces further con¬ 
straints such as the requirement to use a deep lattice po¬ 
tential and a complex band structure even in the weakly 
interacting regime. Consequently, determining optimal 
preparation is beyond the scope of this work. 

To conclude, we briefly outline methods to detect the 
CSL and SFS. Although it is premature to analyze pro¬ 
tocols in detail, as these will depend substantially on the 
specific experimental implementation, it is useful to de¬ 
scribe the basic ingredients that would be required. To 
detect the CSL Ref. [6] suggests methods to probe two 
characteristic properties of topological phases: looking 
for topologically protected, chiral edge currents and in¬ 
troducing a weak attractive optical potential that is lo¬ 
calized to a few lattice sites, which should bind the any- 
onic quasiparticles. Braiding or interfering these quasi¬ 
particles can manifest their anyonic nature. To detect 
the SFS state, one can perform spin-dependent Bragg 
spectroscopy to detect the 2-spinon continuum in the dy¬ 
namic spin structure factor; the most basic signature of 
the exotic nature of this phase is the lack of order and 
existence of gapless excitations. More details of the state 
and its excitations could be revealed by considering more 
structure of the spectrum, similar to that considered in 
Refs. [71]and 72. 
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SUPPLEMENTARY MATERIAL 


Slave rotor mean-field theory: no translational 
symmetry breaking 

Synthetic gauge field case .—Here we give a detailed de¬ 
scription of the slave rotor mean-field theory in the pres¬ 
ence of the synthetic gauge flux. To study the energetics 
as well as the phase transition from the IQH to the CSL 
in Figure [TJb), we first choose the variational ansatz for 
the IQH and the CSL such that tij = £, Jij = J. More¬ 
over, we assume an uniform Lagrange multiplier such 
that hi m h. This is equivalent to replacing the local 


constraint at every site with a global constraint. This 
simplification is justified by the fact that both the IQH 
and the CSL preserve the lattice translation symmetry. 
The rotor and spinon mean-field Hamiltonians are then 
given by 

Hr = - J X + X I- L i+ h + !). (8) 

(Li) i 

Hf = -t X (9) 

(i,j),a i,OL 

where we have replaced the rotor variable e l0i by a uni- 
modular operator such that |<£b| = 1. Since the oper¬ 
ator Li = ^2 a fiif a i — l, h is then thought as a chem¬ 
ical potential. Because (JA L^) = 0, h must vanish for 
the mean-field solutions. With a 27 t/N flux per square 
plaquette and one fermion per site, the Hamiltonian Hf 
gives a spinon band structure with N bands. Only the 
lowest band is filled (for each species) and separated from 
the others by a gap. Moreover, the lowest spinon band 
has Chern number C = 1 for each fermion flavor a. As 
we shown in Table [IJ whether the system is in the IQH 
or the CSL is determined by the behavior of the rotor 
sector. 

To solve the rotor Hamiltonian H r , we implement a co¬ 
herent state path integral formalism in imaginary time. 
We integrate out the conjugate variable Li and obtain 
the partition function that is written as a functional in¬ 
tegration over the variable, 

/’p$tp$p Ae -S-/ T A i (|<I > i | 2 - 1 ). (IQ) 


Here A^ is the Lagrange multiplier introduced to the uni- 
modular condition for the rotor variable at every lattice 
site. The effective action is given by 

S= / 277 X l^r^k| 2 -2J (cos k x +cosky) |$ k | 2 , 

JdT U kGBZ keBZ 

(n) 

where “BZ” refers to the Brioullin zone of the square op¬ 
tical lattice and we have set the lattice constant to unity. 
In a standard spherical approximation for a mean-field 
(or saddle point) analysis, we assume a uniform Lagrange 
multiplier A^ = A. We integrate out the variable 4> and 
obtain the saddle point equation for A, 


1 

A h 


X - = 1 > 

keBZ k 


( 12 ) 


where co>k = [2E/(A — 2J(cosk x + cos ky))] 1 / 2 is the band 
dispersion of the rotor and N s is the number of lattice 
sites. We solve the saddle point equation Eq. (12) self- 
consistently with the spinon mean-field Hamiltonian Hf. 
When the rotor band touches zero energy, the rotor is 
condensed and the internal U(l) gauge field picks up a 



7 


mass due to the Anderson-Higgs mechanism. The ro¬ 
tor and the spinon are then bound together and form a 
fermionic atom. The resulting phase is the IQH. When 
the rotor band is gapped and the rotor is not condensed, 
the internal U(l) gauge field is gapped out by the Chern- 
Simons term and the resulting phase is the CSL. 

In the IQH, the system has N chiral edge modes that 
transport spin quantum numbers as well as atoms. The 
CSL, however, is a Mott insulating state. The atoms 
are localized by the interaction in the CSL. The chiral 
edge states in the CSL only carry spin quantum number 
and cannot transport charge. The effect of the phase 
transition from the IQH to the CSL on the edge states is 
to gap out the mode that transports atoms. 



FIG. 3. The VBS state for the SU(3) model. The spinon 
hopping is zero on the light bonds and is non-zero on the 
dark bond. The tv flux through the bold rectangle is the mean- 
field gauge flux felt by the fermionic spinons in the mean-field 
Hamiltonian H /. The figure is adapted from Ref. 6 

No synthetic gauge field case .—With no gauge flux the 
Hamiltonian is modified by putting <fij —>> 0. For the 
FL and the SFS, we choose the spinon mean-field ansatz 
such that tij = t, Jij = J and hi m h. The spinons 
partially fill the bands and give rise to the spinon Fermi 


surface. Again, whether the system is in the FL or the 
SFS is determined by the rotor sector. Since the FL and 
the SFS only occur for the model without the synthetic 
gauge flux, the rotor sector Hamiltonian is identical to 
Eq. (J8|. When the rotor is condensed, the system falls 
into the FL. When the rotor is gapped, the system is 
in the SFS whose low energy property is described by 
the spinon Fermi surface coupled with a fluctuating U(l) 
gauge field. 

Variational ansatz: VBS states 

As we described in the main text, VBS states become 
favorable in the strongly interacting limit for the SU(3) 
and SU(4) models. As expected, a similar conclusion was 
found in the previous slave-fermion study, i. e. the ground 
state of the SU(3) [SU(4)] Heisenberg model - the U/t 
oo limit of the Hubbard model - on the square lattice 
[Bj favored the VBS state shown in Figure [3] [Figure [ 4 ]. 
In our slave rotor mean-field calculation, we have chosen 
Hf such that the spinons have the same hopping and feel 
the same mean-field gauge fluxes as those ones shown in 
Figure [3] and Figure [4j while the rotor sector has the same 
cluster structure as the spinon sector but experiences a 
different flux (see Table. [T|). 
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FIG. 4. The VBS state for the SU(4) model. The spinon 
hopping is zero on the light bonds and is non-zero on the 
dark bond. The 0 flux through the bold square is the mean- 
field gauge flux felt by the fermionic spinons in the mean-field 
Hamiltonian H /. The figure is adapted from Ref. 0 














